Finding heavy Higgs bosons heavier than 2niw in dileptonic VF-boson decays 
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We reconsider observables for discovering a lieavy Higgs boson (with mu > 2miv) via its di- 
leptonic decays h WW — >■ ivEv. We show that observables generalizing the transverse mass that 
take into account the fact that both of the intermediate W bosons are likely to be on-shell give 
a significant improvement over the variables used in existing searches. We also comment on the 
application of these observables to other decays which proceed via narrow-width intermediates. 



The LHC collaborations have now excluded a heavy 
Higgs boson, with mass between 127 and 600 GeV hi a, 
at the 95% confidence level, assuming that its produc- 
tion and decay proceed according to the predictions of 
the Standard Model. Nevertheless, there is a clear mo- 
tivation to continue the search in this mass region. In- 
deed, theories beyond the Standard Model, with addi- 
tional degrees of freedom, typically predict a lower cross 
section times branching ratio for Higgs production and 
decay in any given channel, for three reasons. Firstly, 
the branching ratio will necessarily be reduced whenever 
there are extra states into which a given Higgs boson may 
decay. Secondly, in theories with an extended scalar sec- 
tor, mixing reduces the cross section to produce a given 
mass eigenstate. Thirdly, since the dominant production 
mechanism for Higgs bosons at the LHC is via loop-level 
gluon fusion, the production cross section will be changed 
by the presence of additional coloured states in the the- 
ory. If one of these states plays the role of cancelling the 
quadratic divergence in the Higgs mass, coming from a 
loop of top quarks, then it necessarily interferes destruc- 
tively with the dominant contribution to gluon fusion in 
the Standard Model, which also comes from a top quark 
loop 0. Furthermore, it is observed (though sadly no 
theorem to the effect has been proven) that the cross- 
section is reduced in all models in which the Higgs boson 
is composite 0-01 • 

One way to search for a heavy Higgs boson is to look for 
evidence of it decaying to the fully leptonic final state Ivliy 
(where I is an electron or muon) which results from an 
intermediate pair of (possibly virtual) W-bosons (sl-fiot. 
A key step in many such analyses comprises a selection, 
or cut, based on a variable that should have a strong 
correlation with the mass of the Higgs boson in 'signal' 
events. Typically this variable will be a transverse mass 
of some kind, of which many have been described in the 
literature [U, ■ The transverse mass variable tot of 
Ref. [13 (see also d, [13) is used in current LHC h -> 
WW — >■ liyli^ searches.^ 



It has previously been observed [l5| that a kinematic 
variable m!^?"'^'^ (which, for reasons that will become 
clear, we refer to henceforth as mij?™'') which demands 
intermediate mass-shell constraints for the tau leptons 
might prove advantageous in /i —>■ rr decays. As was 
indicated in that paper, the same variable (after the 
trivial replacement nir mw) might also prove use- 
ful in the context of ft, — >■ WW when a heavy Higgs bo- 
son has mass greater than 2mw, since the intermediate 
VF-bosons should be produced close to their mass shells. 
Substantiation of this suggestion was left open to fur- 
ther study and one of our intentions here is to perform 
that follow-up. We also intend to go somewhat further, 
in providing a complete discussion of the variables that 
characterize the kinematics of ft — WW — >■ Iflv decays, 
in a sense that will be made explicit below. We will find 
that the pre-existing variables ttit or [lit provide 
a complete description in cases where neither or either, 
respectively, of the intermediate W^-bosons is produced 
on-shcU.^ A complete description of the cases when both 
of the W-bosons are produced on-shell requires not only 
the extant variables m!^'"'"' and mT2 [3, [iSIj but also 
a new variable, which we denote m]^^°'^. Our simula- 
tions suggest, moreover, that 771!^'^^°' gives a significant 
improvement in discovery potential for h — > WW — > li/liy. 

Having finished our preamble, we now describe more 
explicitly how variables of this kind may provide a com- 
plete characterization of the kinematics of events. Since 
the decay of the parent Higgs boson involves invisible 
neutrinos, one cannot find a variable that measures its 
mass directly; rather the best one can do is to find vari- 
ables that bound the mass of the parent in some way. 
Myriad variables of this type may exist and what one 
would like to do, presumably, is to find the variable or 
variables that give the optimal bounds on the mass in an 



^ The variable mx was referred to as mJj?'"" in [Tal to distinguish 
it from other, approximate transverse mass variables in use at 



the time; we now refer to it simply as mx, in line with both the 
ATLAS and CMS collaborations. 

In the mass range < < 2mw, one might reasonably 

assume that one or other of the W bosons is produced on shell 
[TtI . However, simulations reported in |16|| found no significant 
improvement in discovery potential when was substituted for 
rriT in existing ATLAS collaboration searches. 
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event. 

This search for the optimal variables can be done in 
a more-or-less definitive way if one is willing to restrict 
one's attention to kinematics alone. That is to say, sup- 
pose one assumes that a signal event corresponds to some 
decay topology (here the decay of a parent Higgs into 
two intermediate W resonances, followed by decays of 
each intermediate into a combination of visible and in- 
visible daughter particles) . One may then write down the 
various kinematic constraints (corresponding to energy- 
momentum conservation and the mass-shell conditions) 
and ask which values of the a priori unknown masses 
are allowed, in the sense that the kinematic constraints 
then admit solutions in which the unmeasured momenta 
and energies in the event are real and real-positive, re- 
spectively. The resulting allowed mass region is an event 
observable and encodes all of the information which can 
be gleaned about the masses in an event using kine- 
matics alone. The boundary of the allowed region can 
be described in terms of one or more relations between 
the masses, which are themselves event observables and 
which give the optimal bounds on the mass in an event. 
Given multiple events, the allowed mass region is ob- 
tained as the intersection of the allowed regions for each 
event and is itself an observable.'^ 

This abstract recipe leads, in many cases, to simple, 
pre-existing observables. In the canonical example of a 
decay of a parent to a single invisible daughter (together 
with possibly multiple visible particles), the observable 
that results 23 1 is the familiar transverse mass, mx ji^ . 
invoked in the discovery of the VF-boson [l^ in its decay 
to a charged Icpton and a neutrino; for single parent de- 
cays into multiple invisibles, the resulting variable [lij is 
Tot; for identical decays of a pair of parents into invisible 
particles, the variable [23, 26 



271 is TOT2- 



Thus, these variables encode all of the information that 
is available from kinematic considerations alone (subject 
to an assumed decay topology) and there is no point in 
trying to devise further variables to glean further infor- 
mation from kinematics.^ 

In the next Section, we firstly try to prove the equiv- 
alence between the kinematic constraints with one inter- 
mediate on-shell and the variable to^ and secondly the 



^ Whether this observable is the optimal observable for discovering 
the parent is unclear. To settle this, one would first have to 
define what one meant by an optimal discovery observable. For 
discussion, see e.g. |20l422l . 

The equivalence between variables and kinematics is also impor- 
tant for the purposes of determining which particle masses can be 
measured using kinematics alone: as an example, for a (possibly 
pair-produced) parent particle that decays off- shell to multi- 
ple invisible particles, it is a theorem that one can - at least in 
principle - measure the parent mass and the sum of the masses 
of the invisible daughter particles using kinematic information 
alone [MJH. 



equivalence between the constraints with both interme- 
diates on-shell and the variables mlj?^°'' and tot2- Our 
first attempt at a proof succeeds: to^ is equivalent to 
all of the information contained in the kinematic con- 
straints when one intermediate is on-shell. But for the 
topology with both intermediates on shell, we fail. In- 
deed, we show that toJJ?^'''' and tot2 do not capture all 
of the information in the kinematic constraints. Rather 
there is a third, distinct variable, which we call to!^^^°'^, 
and which gives an upper bound on the mass of the Higgs 
(or whatever parent particle is being considered). 



In what follows we show that while 



T lower 



only 



marginally enhances the discovery potential for a Higgs 
boson using current ATLAS collaboration search strate- 
gies, m^^'^'^ gives a significant improvement. In an 
Appendix, we discuss issues related to the existence of 



to!S™'' and ml^P"'. 



KINEMATIC EQUIVALENCE 

Since we are concerned with the three transverse 
masses tot, rnl^ and to!^™'', it is helpful to recall and 
compare their definitions. All three make the assump- 
tion that the observed Ivlv final state resulted from the 
decay of a single progenitor particle or resonance (i.e. the 
Higgs in the case of the signal) whose mass is, a priori, 
unknown. All three assert that the sum of the trans- 
verse momenta of the neutrinos should equal the observed 
missing transverse momentum. All three accept that: (i) 
the way the transverse momentum is shared between the 
two neutrinos is unknown, and (ii) that the longitudinal 
momentum component of each neutrino is unknown. Fi- 
nally, all three transverse masses are defined to be the 
greatest possible lower bound for the unknown mass of 
that progenitor particle, subject to consistency and any 
remaining constraints.^ It is only in these remaining con- 
straints that the transverse masses differ: tot applies 
none, to^ permits only neutrino momenta that place at 
least one VF-boson on mass-shell, while TOlj?™°'' allows only 
neutrino momenta that place both VF-bosons on mass- 
shell. In short, history would have been kinder to us 



had Tot, rn1^ and to!^^*"" been named 
and respectively. The variables and constraints 

are summarized in Table [D 

We now attempt to prove the equivalence between the 



Equivalcntly, they may be defined to be the smallest possible pro- 
genitor mass obtainable by a search over all unknown neutrino 
momentum components that satisfy the remaining constraints. 
Indeed, this brute force search technique is used in the authors' 
implementations of and mij?™""^ calculators. These C++ im- 
plementations, free of external dependencies, are available on re- 
quest. A simple analytic answer is available for jtit - sec Eq. (f ) 

of fii. 
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Name Ref. Bounding Constraints 

mT \8, 13, 14} H^'Hf, © to Q 

rriT [16] Hf'H^ O to © and either (g]) or © 

flower JJM^^ Q tO © 



(a) Variables (b) Constraints 

TABLE I: A comparison of the three mass-bound variables. The final state four-momenta of the two visible leptons are p'^ 
and respectively. The hypothesized four- momenta of the two neutrinos are given by the two p'^. The missing transverse 
momentum is p^. In all three cases, the variable named in (a) is defined to be the maximal lower bound on the invariant 
H'^H^, where H'^ — p'i +pf + p'i^ '^p't ^ subject to the indicated subset of the constraints shown in (b). 



p,.p^,^ = PiPifj. 


= 0, 


(1) 


p,!rp„^ = Pi Pit, 


= 0, 


(2) 


PuT + p'vT 


= fj,, 


(3) 


{Pu +Pt){Pi^ti +Plp.) 


2 


(4) 


{pT+pT){p'^+pi,) 


2 

= m^. 


(5) 



variables and mlp™*"' and the corresponding kinematic 
constraints. 

One intermediate on-shell 

The theorem to be proved is a formalization of the 
colloquial: to^ is equivalent to all the information con- 
tained in the kinematic constraints. To wit, in a notation 
in which = {E,pT,q)- 

Theorem 1 The kinematic constraints ^ to ^ and 
either ^ or ^ admit a solution with momenta in K 
and energies in M"*" ijf. 

rufi > m^. (6) 

Here, is defined as the minimum value of 

rrih > (7) 

obtained by varying the unobserved momenta in R and 
energies in R+. This definition makes the necessity of 
ruh > nfi^ obvious; it is the sufficiency of the condition 
that requires deliberation. 

So, does ruh > imply that the kinematic con- 
straints have a solution? One way to show this would 
be to consider every ruh > mt^ and explicitly construct 
a solution therefor. An easier way is to show that one 
can find solutions corresponding to arbitrarily large ruh 
and then to invoke continuity and the intermediate value 
theorem to show that one can find solutions for any m/j 
between and oo. 

Let us then exhibit a solution for arbitrarily large m/j 
that satisfies constraints ([T]) to (U) (if the other W is 
on-shell, the required solution can be obtained by inter- 
changing primed and unprimcd quantities). It is given 
by 

= YimJ^ql^+fl^f^, ql), (8) 
P'l = i\q.l0,q.) (9) 



where 5^ is given by either of 

= ^ ± ^if+Pfr) • (10) 
This yields a value of given by 
ml = Hm (m^ 4- 2{pi + p^) ■ p'l + 2{p[ + pi + Pu) ■ 

— >oo 

-^2{{E[ + Ei + \q,\)^ql^+fl 

- (pit+pit) ■ fj, - {q'l + qi + qy)q'v)) 

2{{E[ - q'l) + {El - qi) + {\q,\ - q,))q'„ (11) 

where in the last line we assumed > without loss 
of generality (this amounts to a choice of which proton 
beam constitutes the +z direction). The terms in paren- 
theses in the coefficient are each positive semi-definite 
and so no cancellation of the coefficient of q'^ can occur. 
Thus ruh grows without bound unless all terms vanish. 
To do so, both leptons must have vanishing transverse 
momenta, in which case they would fall outside the de- 
tector acceptance and this would not be identified as a 
di-lepton plus missing momentum event. | 

Both intermediates on-shell 

It is now easy to show that an analogous theorem can- 
not be proven for mlj?"°'', which applies to the topology 
in which both constraints Q and ^ are satisfied. The 
reason is that one cannot find solutions to the kinematic 
constraints corresponding to arbitrarily large ruh- In- 
deed, one can easily convince oneself that to get arbitrar- 
ily large rnn, at least one of the neutrino momenta must 
become arbitrarily large. But each neutrino is now pro- 
duced in the decay of a W-boson, which simultaneously 
results in a charged lepton of finite (and measured) mo- 
mentum. Now, there is only one way in which a two-body 
decay can produced one daughter of finite momentum 
and another with infinite momentum in the lab frame: 
the rest frame of the W must be infinitely boosted exactly 
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in the direction of motion of the neutrino (as measured 
in the W rest frame). This doesn't work, because the 
lepton momenta are then not only finite but arbitrarily 
small, in contrast with their measured values in a generic 
event. 

Thus, there is more information in the kinematic con- 
straints than is captured by ml^^'^'' alone. In particular, 
since cannot reach arbitrarily large values in events, 
there is a distinct variable, obtained by maximizing the 
sum of lepton and neutrino momenta in events, subject 
to the above constraints. This variable will be bounded 
below by to/j and we call it m^^°'^ . By the intermediate 
value theorem, this variable, together with mJj?™'' docs 
contain all of the information in kinematics. 



MONTE CARLO SIMULATIONS 

Simulations, similar to those performed in Ref [l3| . 
are used to test the extent to which using mJj?™°'' and 
n%^^'^'^ might be expected to enhance the statistical sig- 
nificance of a /i — WW signal above the dominant WW 
continuum background and the ti background. We use 
the HERWIG 6.505 [H ^ Monte Carlo generator, with 
LHC beam conditions {^/s = 7 TcV). Our version of the 
generator includes the fix to the h — > W^W^'-*'' spin corre- 
lations described in [sJl ■ We generate unweighted events 
for Standard Model Higgs boson production (gg h) 
and for the dominant backgrounds, qq — > WW, and tt 
production.^ We use the leading order Standard Model 
cross section for all values of nih- 

The missing transverse momentum is calculated from 
the negative sum of the px of visible particles within the 
fiducial region pr > 0.5 GcV and |?7| < 5. The detec- 
tor resolution is simulated by smearing the magnitude of 
the missing momentum vector with a Gaussian resolu- 
tion function of width ap /f^ = 0.4 GeV^/^/-\/S where 
S is the sum of the |pt| of all visible fiducial particles. 

For each value of , fifty pseudo-experiments are per- 
formed, each corresponding to 10fb~^. Selection cuts are 
applied, requiring: 

• Exactly two Icptons t G {e, /x} with > 15 GcV 
and |?7| < 2.5 

• Missing transverse momentum, > 30 GcV 

• 12 GcV < mu < 300 GeV 

• No jet with pr > 20 GeV 



m Z TT rejection: the event was rejected if \mT-T ~ 
mz\ < 25 GcV and < < 1 for both i e {1, 2}^ 

For each of the three candidate discriminant variables 
mT, mJj?™"^, and mJJP^"'^, and for each trial value of nih, a 
likelihood is calculated for each of two hypotheses H: a 
background-only (B) hypothesis and a signal plus back- 
ground hypothesis {S + B). The log likelihood 




is calculated from the average over 50 pseudo- 
experiments. For each pseudo-experiment j the num- 

(1) 

ber of events observed in the ith bin is n~ . The ex- 
pected number of events x]^ for each bin and for each 
of the two hypotheses H is calculated from indepen- 
dent Monte Carlo samples containing 10^ events for each 
trial value of mh and 10^ events for the WW back- 
ground. The likelihood for each np-* is calculated us- 
ing a Poisson distribution with mean . The difference 
— 2Alog£ = —2 (log£s+B — log^Cs) quantifies the sta- 
tistical significance with which we can expect the Higgs 
boson signal to be observed. 

The discovery potential using separately tot, m}^™°'^, 
or ?7i!Jf is shown in Figure [T] The first result is 
that when rrih ~ 2mw the doubly-on-shell lower bound 
variable, mJ^?™"^ has significantly better sensitivity than 
the simple transverse mass tot (or indeed the singly- 
on-shell variable m^, which is almost indistinguishable 
from Tot [lij). Secondly, we observe that the distribu- 
tion of the doubly-constrained upper bound m^^'^^ shows 
a markedly greater discovery potential again, relative to 
either tot or mJj?™°'^. Since any event satisfying the con- 
straints of TableUwill generate both a lower and an upper 
bound, and since those bounds are not closely correlated 
(Figure [2]) each adds information relative to the other. 

OTHER FINAL STATES 

Wc note that it is also possible to construct equivalent 
bounds for any other decays of the form A ^ B+C where 
B and C are produced near their mass shell and each 
decays to an invisible particle and one or more visible 
particles. A case relevant to LHC Higgs physics is the 
decay mode h — > rr, followed by semi-invisible decay of 
each of the r leptons. Previous simulations [lH indicated 
that adding the tau mass constraints led to improved 
discrimination of the signal from the dominant Z — > rr 



Other backgrounds, such as Z — tt, are rendered sub-dominant 
by the cuts discussed below. 



The variable Xi is the momentum fraction of the ith tau carried 
by its daughter lepton and rriT-r is the di-tau invariant mass. 
They are calculated using the approximation that each t was 
coUincar with its daughter lepton. 
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FIG. 1: Higgs boson discovery potential as a function of nih, 
for each of three different kinematical variables, ruT, mx"""^, 
and 771^'^'^'^. The centre of each band indicates the differ- 
ence in log likelihood between models with and without a 
Higgs boson contribution. Lower values correspond to better 
discovery potential. The half-width of each band gives the 
root-mean-squared variation over the 50 pseudo-experiments. 
The integrated luminosity simulated is 10 fb"'^. 



CONCLUSION 

When a heavy higgs boson with mass > 2mw decays 
via h — )• WW — > £u£i/, we have shown that one can use 
knowledge about the narrow width of the intermediate W 
bosons to construct both upper and lower bounds on mh- 
We demonstrate that these two bounds together make 
maximal use of the kinematic information about ruh in 
each event. Since the two bounds are reasonably uncorre- 
lated, the ATLAS and CMS experimental collaborations 
should be able to obtain improved kinematic sensitivity 
in an analysis which simultaneously makes use of both 
bounds. 

We note that the same suite of variables can find use 
for discovery or mass determination in other decays of 
the form A ^ B + C followed by semi-invisible decay of 
each narrow-width intermediate particle. 

This work was supported by the Science and Technol- 
ogy Facilities Council of the United Kingdom, by the 
Royal Society, by the Institute for Particle Physics Phe- 
nomenology, by Merton College, Oxford, and by Peter- 
house, Cambridge. We gratefully acknowledge T.J. Khoo 
for helpful comments. 



Appendix: Remarks on the existence of m.^^"^ and 




150 



160 170 180 



190 200 210 220 
flower / Q^y 



FIG. 2; The distribution of m^'''""' 
boson signal sample with ~ 200 GeV. 



m!r°"'="' for a Higgs 



background, using the lower bound ml^™'' 15[. However 
further simulations indicate that there seems to be little 
additional discrimination available from the upper bound 
in this case, perhaps because of the large hierarchy in 
masses between rrih and mr- 



Note that there exist 'input momenta' for which it is 
not possible to satisfy all the constraints required by 
j.^iower 77i!^PP°''j to which the following arguments 

equally apply). For events containing such momenta, a 
value of mi^'*'°'^ cannot be computed. In effect, the vari- 
able does not exist as those events are incompatible with 
the assumptions at face value. ^ It is known [l^ that 
m}^'"''' exists if and only if mT2 [H, for the dilepton 
system is less than mw ■ This condition may easily be vio- 
lated if one or both Il^-bosons is off mass-shell, confirming 
that there is little sense in applying m}^'"°'^ to events with 
jTift < 2mw- The criterion mT2 < can also be vio- 
lated in mh > 2mw events with sufficiently large detector 
mis-measurements, even if both M^-bosons are produced 
on-shell. The largest source of mis-measurement is the 
missing transverse momentum, whose distribution can 
have a large, non-Gaussian tail. Though not discussed in 
[iB l. this potential problem might be thought more wor- 
risome in the case of /i rr decays, since the typical size 
of an LHC mismeasurement greatly exceeds m^.. It is 
thus hard to see, naively, how the condition 771^2 < 
may be maintained. Nonetheless, one is assisted by a 



Note also that the variables mx and exist for all input mo- 
menta, and are therefore usable even when m}^"^'^^ is not. See 
diseussion in llSll . 
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property of the mT2 variable under which large mis- 
measurements of can lead, in a non-negligible number 
of events, to small changes (relative to rrir) in the value 
of mT2i provided that the intermediate particles in the 
production process are sufficiently light. Indeed, in the 
limit that the intermediate pair of particles are massless, 
their decay products are collinear and so the true vec- 
tor lies in the smaller of the two sectors bounded by the 
two visible transverse momenta. It follows that for such 
configurations mT2 must vanish, as there exists a split- 
ting of the (in fact the 'true' splitting) for which the 
transverse mass of each pair of decays is zero. That tot2 
is 'small' for well-measured events in this limit is not sur- 
prising - after all, mT2 is supposed to be bounded above 
by the mass of either member of the pair, and we are 
taking the limit in which those masses go to zero. What 
is surprising, if anything, is that 771^2 can be forced ex- 
actly to zero in a large number of events, even for finite 
yet small intermediate particle masses. In this sense we 
can say that mT2 goes to zero faster than mr- It fol- 
lows that for such signal events, any mis-measurement 
of (even large ones) that keeps it 'between' the visi- 
ble particles' transverse momenta, loads to no observable 
change in mT2- Moving away from the above limit, such 
that the intermediate particles now have non-zero but 
still small intermediate masses, we see that changes in 
771T2 can still be small, even for large mis-measurements 
of , provided that mis- measured remains in the ap- 
propriate region of the transverse plane. 

This goes some way towards explaining why it was 



found in [l5| that rn}™" could remain well-defined for 



a large number of h ^ tt events, even in the presence 
of momentum mis-measurements. For h — !■ WW events 
this issue is much less of a concern; we find that mij?'"'''' 
exists for all but a few per cent of the h — >■ WW signal 
events we simulate. Wc attribute the greater frequency 
with which mij?™'' exists to the greater ease with which 
the momentum mis-measurement errors (which have a 
natural scale far below mw) may be incorporated into 
the intermediate particle mass constraints (whose natu- 
ral scale is mw)- 
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